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Abstract 

As the new results for the massive Thirring model the L-matrix and the al- 
gebraic relations for its action angle variables are given. So it is shown most 
directly that this model which describes self-interacting relativistic Fermions in 
one-dimensional space is a quantum integrable system. 
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1 Introduction 



In this note we present the L-matrix for the massive Thirring model [1] defined 
by the Hamiltonian 

pM 

H= dx{4> + (-iaz)d x (i) + m(i) + ai4>-2siiy{2c)(i)'t4>t4>24>i}, (1) 
Jo 

where Oi are the Pauli matrices. The Fermion fields 4>{x) = (0i(x), 02 (x)) satisfy 
the usual equal time anti-commutation relation {0 M (a;), (f>t(y)} = S^, u 5(x — y), 
the indices refer to right or left moving particles. The constant c < ^ gives the 
strength of interaction between the Fermions. As a solvable and simple model the 
massive Thirring model has attracted much interest (see refs. [2-11]). 

The L-matrix corresponding originally to the Lax-pair for nonlinear partial 
differential equations is an operator that determines a " canonical" transformation 
of the field variables (</> + , 0) to variables of action-angle type, which are the matrix 
elements of the so-called monodromy matrix. The integrability can be shown by 
constructing a i?-matrix, which gives the information about the algebraic relations 
between these action-angle variables. Hence the L-matrix and the i?-matrix are 
crucial for the study of an integrable system. 

The Thirring model is equivalent to the quantum sine-Gordon model [2]. For 
this bosonic model the L-matrix and i?-matrix have been given by Sklyanin, 
Takhtadzhyan and Faddeev [4] . We think it is an interesting problem to construct 
these matrices also for its fermionic counterpart especially since supermatrices are 
the tools which are not so common with integrable systems. 

Here we will study the L- and L'-matrices for the Thirring model. The method 
we used originally is to find the L-matrix by taking the continuum limit from the 
.R-matrix of an inhomogeneous six vertex model. It is known that the six vertex 
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model can be used to construct the Bethe ansatz solution of the massive Thirring 
model [8] and that the Hamiltonian should be connected to a six vertex model 
with staggered weights [10]. Actually, we used this method before to derive the L- 
matrices for the sine-Gordon model [13] and a bosonic system with the hamiltonian 
similar to equ. (1) [12]. Since in the present case the derivation is lengthy and 
similar to the one in [12], we use here a different approach. Given the result, we 
have to prove that we have found the L-matrix of the Thirring model. So we have 
to find the .R-matrix starting from the L-matrix in order to show the integrability. 
Since L-matrix determines the monodromy matrix, we have to check whether the 
latter contains really the Hamiltonian given above as one of the simple conserved 
entities. 



2. L-matrix 

The action-angle variables or more specific the monodromy matrix T(u) = 
T(x = M\u), M is the length of the system, are defined by the differential equation 

dT(x\u) r . 

K -^- =: L(x\u)T(x\u) : (2) 

clx 

and the boundary condition T(x = 0\u) = I, the identity matrix. The colons mean 
normal ordering of the Fermi operators and the u is the spectral parameter. 

We found for the L-matrix after taking the continuum limit and introducing 
fermion fields 

L(x\u) = i-msinh-uT3 + + S(x\u). (3) 

Here we distinguish the Pauli matrices in classical space r\ and those in quantum 

space Oi used before. For E and S we have 

(-<P + (x)(l-e^Mx), \ 

[X) V 0, -d>+(x)(l + e-^)d>(x))' [) 



S(x\u) = ivmsinc _«, + / \ , «, + , x n • (5) 

v 1 ' e 2 0^(x) + e 2 / 

We have to choose the L-matrix as a super-matrix with its row or column parities 
= and p(2) = 1, consequently T(x|u) as a solution of equ.(2) is also a 
super-matrix. 

The integrability of the Thirring model can be shown by finding the Yang- 
Baxter equation for the monodromy matrix T(u) 

R(u - v)T(u)sT(v) = T(vfsT(u)R(u - v), (6) 



where the tensor product indicated by the symbol s is of the super- form (AsB) i '• 
= (-l)pO')(p(*)+p(0)^ ijfcj B J . i ,. The matrix R has the same form as for the six vertex 
model, note however that the spectral parameter u is imaginary in the usual six 
vertex model 



R(u) = 



/sinh(f + ic), 0, 

0, sinhic, sinh-|, 

0, sinh^, sinhic, 

V 0, 0, 0, sinh(f+ic) 



(7) 



In the following let us prove the Yang-Baxter equation (6). We can rewrite the 
differential equation (2) as an integral equation for < x < M 



T(x\u) =e ^sinW 3 

+ / dze^ {x - z)sinhuT3 : (£(z) + S(z\u))T(z\u) :, 
Jo 

from which we can get more easily two auxiliary equations 

<j) v (x)T{x\u) = t 3 T{x\u)t 3 4> v {x) + ^ : E v (x\u)T(x\u) :, 

I 

T(x\u)(/>+(x) = 4>+(x)t 3 T(x\u)t 3 + - : F u {x\u)t 3 T{x\u)t 3 :, 



(8) 



(9) 



with 

E v {x\u) = -<l> v {x)(\ - r 3 e iCT3( - 1)1/+1 ) + 'e { ~ iy r~ 

(10) 

F v (x\u) = -<f>t(x)(l - r 3 e lCT3( - ir+1 ) - iV^^(-l) v e^" %r+. 

Defining the tensor product in the Yang-Baxter equation (6) as 

K(x\u, v) = T(x\u)sT(x\v), (11) 

one can get using the last equs. (9) and (10) a differential equation for this tensor 
product 

d x K(x\u,v) =: L(x\u,v)K(x\u,v) : (12) 
with a L-matrix depending now two spectral parameters u and v 

L(x\u, v) = L(x|u)f 1 + l'sL(x\v) + F v (x\u)% E v (x\v). (13) 

v 

One can show only by an explicit calculation that 

R(u — v)L(x\u, v) = L(x\v, u)R(u — v). (14) 

The equation means that i?-matrix can exchange u and v in L(x\u,v) in manner 
the Yang-Baxter relation postulates. The equs. (11-13) above defining the tensor 
product K(x\u, v) from L(x\u, v) are of course also valid if the spectral parameter 
u and v are exchanged. This implies that an equation like equ. (14) must hold 
also for K 

R(u — v )K(x\u, v) = K(x\v, u)R(u — v), (15) 
which is the Yang-Baxter relation (6) taking x = M. 

Hence the transfer matrix t(u) = T(u)u — T(u) 2 2 which is the super-trace of 
the T(u) must commute for different spectral parameters, i.e. 

[t(u), t(v)] = 0. (16) 
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So the definition of the L-matrix (3) generates a quantum integrable system. The 
main problem is now to show that this quantum system is the massive Thirring 
model. 

3. Hamiltonian and momentum 

Here we will find the Hamiltonian of the Thirring model (1) from the transfer 
matrix t(u). This shows directly that the L-matrix (3) gives the Thirring model. 
Using an integral form of equ. (2) similar to (8) , however expanding T with respect 
to S given by equ. (5), one has 

f(x\u)=:Q:+ [ dze- imzT3Sinhu : Q(z)S(z\u)f(z\u) :, (17) 
Jo 

where f(x\u) = e ~ l ^ XT3 sinhu T(x\u), Q(z) = exp(f* dyZ(y)) and Q = Q(z = 0). 
In order that such an expansion makes sense one must add an imaginary part to 
the spectral parameter, for example, u — > itizf to have simple expressions. Thus 
from 

T(x\u ± i-) =e i ^ XT3 sinh ( u±i f ) :Q:+... 

2 ' (lg) 

— e Tf IT 3 coshu . q . _|_ 

we can see that the transfer matrix t(u ± if) decrease or increase rapidly for 
u — > ±oo. After iterating equ. (17) we see that the expansions 

f(x\u-i-) u =: Qn : +msinc / dz x \ dz 2 e- m ^ Zl - z ^ coahu 

2 " Jo Jo ~ (19a) 

x:<f>_(z 1 )<f>±(z 2 )Q 11 :+--- 

and 

f(x\u + i?-) 22 =: Q 22 :+msmc f dz x C dz 2 e~ m ^-^ coshu 

2 " J J Q - (196) 

x : ^+( Zl )^ + (z 2 )Q 22 :+■■■ 

have nontrivial limits for u — > ±oo, where 

^±{z) = Q{z) 11 {e-^4> 1 {z)-e^- ( t )2 {z))Q-\z) 22 



and 

^+(z)=Q(z) 22 (e-^<P+(z)-e^<P+(z))Q- 1 (z) 11 . 

By choosing the sign of the imaginary part of u we have picked out the exponential 
growing contributions. Neglecting the decreasing part related to exp(— e ±u ) for 
u — > ±00, the transfer matrix multiplied by exp{±i^Msinh(« ± if)} is simply 
T(u — i\ ) 11 or — T(u + 1^)22 in this limiting case. We combine the contributions 
defining the generators 

G±(u)= lim {e ±ic \n[e- i ^ Msinh ^- i ^t(u-il) 

-e^ c ln(e^ Msinh{u+ ^h(u + i-)]}. (20) 

The generators G± can be calculated from equs. (19) using partial integrations 
repeatedly. In this way one obtains a series in e^ u for u — > ±00 

G±{u) = Y,C± s e* su . (21) 

s>0 

The calculation is tedious but straightforward. Every term of the expansions Tn 
and T 22 has contribution to G±, even to the first order C±\. Fortunately, all 
coefficiences of the factor e ±u can be summed up and give us the wanted results 

8 sin c Tn 

C +1 = / dz[i</>td z <h + Trtttfa + 02 0i) - sin(2c)^40 2 0i)] (22) 

m J - 2 

and 

s sin p z*-^ 

C-l = / dzi-itfdzfa + "(0^2 + 0^1) - Sm(2c)<P+ 02^0201)]- (23) 

m Jo 2 " 

It can be seen that the generators G_ and G + give the conserved quantities for 
right and left moving particles respectively. The Thirring system includes both 
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contributions of right and left particles. The sum of the two first order coefficients 
is just the Hamiltonian of the Thirring model (1) whereas the difference is the 
momentum 

pM 

P = -i dz(p + (z)d z (j)(z). (24) 
Jo 

The other coefficients C± s (s = 0, 2, 3, • • •), if they are not zero, should give other 
conserved quantities of the Thirring model. These are not easy to calculate and 
cannot be studied here. 

4. Conclusion 

In this note we have described the L- and /^-matrices for the Thirring model 
which were not given before for a fermionic relativistic theory. For quantum inverse 
scattering transformations the L- and i?-matrices are important operators, for 
example, for studying the algebraic Bethe ansatz and the inverse problem of the 
Thirring model in the sense of the works [4] and [14]. Also it is interesting to 
note that the Yang-Baxter equation (6) for the monodromy matrix has a super 
structure, i.e. the diagonal elements of the monodromy matrix are of bosonic type, 
whereas the off-diagonal ones are of fermionic type. Hence equation (6) gives us 
also a graded and deformed algebra or a graded quantum algebra (see refs. [15,16]). 
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